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ABSTRACT 

Gravitational waves (GW) generated during a core-collapse supernova open a 
window into the heart of the explosion. At core bounce, progenitors with rapid core 
rotation rates exhibit a characteristic GW signal which can be used to constrain the 
properties of the core of the progenitor star. We investigate the dynamics of rapidly 
rotating core collapse, focusing on hydrodynamic waves generated by the core bounce 
and the GW spectrum they produce. The centrifugal distortion of the rapidly rotating 
proto-neutron star (PNS) leads to the generation of axisymmetric quadrupolar oscil¬ 
lations within the PNS and surrounding envelope. Using linear perturbation theory, 
we estimate the frequencies, amplitudes, damping times, and GW spectra of the os¬ 
cillations. Our analysis provides a qualitative explanation for several features of the 
GW spectrum and shows reasonable agreement with nonlinear hydrodynamic simu¬ 
lations, although a few discrepancies due to non-linear/rotational effects are evident. 
The dominant early postbounce GW signal is produced by the fundamental quadrupo¬ 
lar oscillation mode of the PNS, at a frequency 0.70 kHz < / < 0.80 kHz, whose energy 
is largely trapped within the PNS and leaks out on a ~ 10 ms timescale. Quasi-radial 
oscillations are not trapped within the PNS and quickly propagate outwards until 
they steepen into shocks. Both the PNS structure and Coriolis/centrifugal forces have 
a strong impact on the GW spectrum, and a detection of the GW signal can therefore 
be used to constrain progenitor properties. 
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1 INTRODUCTION 


Rotating iron core collapse in a massive star (M > 8Mq, 
resulting in a core-collapse supernova [CC SN]) was one of 
the first potential sources of gravitational waves (GWs) con¬ 
sidered in the literature ||Weber||l966[ |Ruffini fc Wheeler| 


|1971 see Ott 2009 for a historial overview). GWs are of 


lowest-order quadrupole waves and rotation naturally drives 
quadrupole deformation (oblateness) of the homologously 
{v oc r) collapsing inner core of a rotating massive star. 
When the inner core reaches nuclear densities, the nuclear 
equation of state stiffens, stopping the collapse of the in¬ 
ner core. The latter overshoots its new equilibrium, bounces 
back (a process called “core bounce”) and launches the hy¬ 
drodynamic supernova shock at its interface with the still 
collapsing outer core. Subsequently, the inner core rings 
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down, shedding its remaining kinetic energy in a few pul¬ 
sations, then settles to its new postbounce equilibrium and 
becomes the core of the newly formed proto-neutron star 
(PNS). The entire rotating bounce-ring down process in¬ 
volves rapid changes of the inner core’s quadrupole mo¬ 
ment and tremendous accelerations. The resulting GW burst 
signal has been investigated extensively both with ellip¬ 


soidal models (e.g., Saenz &i Shapiro 19781 and with detailed 


multi-dimensional numerical simulations (e.g., |Muller|1982 

Monchmeyer et af.||19911 Yamada & Sato||1995 |Zwerger &| 

Mufler||1997| Uimmefmeier et af. 2002 

Kotake et al. |20031 

Ott et al. 

2004 Dimmelmeier et al.| 

20081 Obergaulingerl 

et al.|2006 

Ott et al.|2007 2012 1 . 


Based on this extensive volume of work, it is now clear 
that rotating core collapse proceeds mostly axisymmetri- 
cally and that nonaxisymmetric dynamics sets in only within 
tens of milliseconds afte r bounce (|Ott et al!]|2007| |Schei-| 
degger et al.||2008 2010 Kuroda et al.||20i4 |. Only rapidly 
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2 J. Fuller et al. 


rotating iron cores (producing PNSs with central spin pe¬ 
riods < 5 ms) generate sufficiently strong GW signals from 
core bounce to be detected throughout the Milky Way by 
Advanced-LIGO-class GW observatories (Aasi et al. (LIGO 
Scientific Gollaboration) 2014 Ott et al. 2012 hereafter 
012). Since the cores of most massive stars are believed to 
be slowly rotating at core collapse (> 90%, e.g., Heger et al. 
2005[[0tt et al.|2006 Langer|2012 ), the detection of GWs of 
a rotating core collapse event may be exceedingly rare and 
GW emission in CG SNe may be dominated by neutrino- 


driven convection instead (e.g.. 

Muller & Janka|1997| 

Muller 

et al. 112004 

|Ott|2 

0091 |Kotake||20131 |Ott et al.||20131 

Muller 

et al.|2012| 

2013] 

Murphy et al.|2009|. However, if a 

rapidly 

could 

rotating core collapse event were to be detected, i 


possibly be finked to an energetic CC SN driven by mag 
netorotational coupling (e.g., Burrows et al.|2007 Takiwaki 


et al. 2012| Mosta et al. 20141. Abdikamalov et al. (20141 


(A14 hereafter) have shown that the angular momentum of 
the inner core can be measured from the observed GW sig¬ 
nal. 

The morphology of the GW signal from rotating core 
collapse, bounce, and ring down is uniform across the en¬ 


tire parameter space of plausible initial conditions (Dim- 
melmeier et al. 2007 2008[ 012) and essentially independent 
of progenitor star mass. It consists of a first prominent peak 
associated with core bounce (the “bounce signal”, cf. Fig- 
ure[^ and a lower-amplitude but longer duration oscillatory 
ring down signal, persisting for ~ 20 ms after bounce. The 
ring down signal is peaked at a GW frequency of 0.7—0.8 kHz 
(somewhat dependent on equation of state and rotation rate, 
Dimmelmeier et al.|2008| 012; ?, in prep, hereafter K15) and 
may be correlated with variations in the early postbounce 
neutrino luminosity (012), suggesting that the ring-down 
oscillations of the PNS (here defined as the inner 20 km of 
the postbounce star) are connected to the excitation of a 
global PNS oscillation mode at core bounce (012). 

While the GW signal of rapidly rotating CG SNe can be 
straightforwardly computed from complex nonlinear multi¬ 
dimensional hydrodynamic simulations, its features are not 
understood at a fundamental level. The goal of the present 
investigation is to provide such a basic understanding of the 
signal features. To do this, we employ semi-analytical, lin¬ 
ear calculations of the wave-like fluid response produced by 
the bounce of the inner core. These calculations shed light 
on the physical mechanisms responsible for the GW signals 
discussed above, and their simplicity complements the com¬ 
plexity of the simulation results. However, our methods only 
provide a qualitative explanation for GW signals, the sim¬ 
ulations are needed for precise quantitative predictions. Al¬ 
though oscillations of PNSs have previously been examined 
in several papers (see, e.g., Ferrari et al.||2003 20041, these 
works focus on PNS oscillation modes well after (> 100 ms) 
bounce, and they do not investigate the physics of waves 
excited by the bounce itself. 

We find that the postbounce fluid response has a few 
distinguishing characteristics. First, the bounce excites a 
train of radial, outwardly propagating acoustic waves. Be¬ 
cause the background structure is centrifugally distorted by 
rotation, these waves are only quasi-radial and contain a 
quadrupole moment, allowing them to emit the GW respon¬ 
sible for the bounce signal discussed above. Second, the cen¬ 
trifugal distortion of the progenitor leads to the excitation 



Figure 1. Top: Central density pc of the rapidly rotating A3O05 
simulation of A14. Time is measured from the moment of bounce, 
defined as the moment when pc peaks. Bottom: GW strain from 
the same simulation, computed at a distance D = lOkpc. The 
large amplitude signal near bounce is the “bounce signal” cre¬ 
ated primarily by the quasi-radial bounce of the centrifugally dis¬ 
torted inner core (012, A14). The smaller amplitude “ring down” 
oscillations after bounce are created primarily by quadrupolar os¬ 
cillations of the PNS. 


of a train of axisymmetric quadrupolar waves. Some of these 
waves are reflected at the edge of the PNS, causing them to 
interfere to create standing waves, whose energy is peaked 
near the oscillation “modes” of the PNS. The dominant GW 
signal is produced by the fundamental PNS oscillation mode 
at a frequency /f_mode ~ 0.75 kHz, accounting for the most 
prominent peak in the GW ring down signal. The GW signal 
damps on 10 ms timescales as the wave energy leaks out of 
the PNS into the surrounding envelope. 

Our paper is organized as follows. In Section we in¬ 
troduce our semi-analytical framework for calculating the 
spectrum of waves excited at core bounce, and we discuss 
the properties of the resulting waves. Section investigates 
the subsequent wave damping, and discusses the complica¬ 
tions introduced rotational and relativistic effects. In Sec¬ 
tion 1^ we present the GW spectra produced by the waves 
and compare with the GW spectra seen in simulations. We 
conclude in Section with a discussion of our results and 
future avenues for theoretical development. 
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Supernova Seismology 3 


2 OSCILLATIONS EXCITED AT BOUNCE 


As described above, the GW spectrum of a rapidly rotat¬ 
ing supernova near core bounce consists of a bounce signal 
and a ring down signal. The bounce signal has a short dura¬ 
tion in time and is thus characterized by a broad frequency 
spectrum, while the ring down signal lasts longer and has a 
spectrum peaked at discrete frequencies. Our main goal is 
to understand the physics of the ring down signal, although 
our methods also shed some light on the spectrum of the 
bounce signal. 

During CC, the inner core of a massive star progenitor 
collapses into a PNS, while the outer core forms the shocked 
region surrounding the PNS core during the on-going su¬ 
pernova. In this paper, we refer to the inner r < 30 km 
{p ^ 10^^gcm“®) as the PNS, while the low density sur¬ 
rounding regions are the envelope. We choose this definition 
because the bounce excited waves can become trapped in the 
inner ~30km. However, note that at a few ms after bounce, 
the PNS only has a mass of Mpns ^O.GMq while the mass 
within the inner 300 km is M ~ IMq, and we therefore ex¬ 
pect all of the material within our computational domain to 
eventually accrete onto the central compact object. In the 
brief (less than a second) period following the PNS bounce, 
but preceding the supernova, the inner ~100 km (i.e., re¬ 
gions below the shock radius) of the supernova core is in 
approximate hydrostatic equilibrium ( Janka|2001 1. 


2.1 Models 


To generate background models for wave excitation and 
propagation within the postbounce supernova structure, we 
use simulation outputs generated by A14. These simulations 
are run with the 2-dimensional version of the CoCoNuT code 


I Dimmelmeier et al. 2002 2005) in axisymmetry and con¬ 


formally flat general relativity. These approximations are all 
appropriate for understanding the axisymmetric waves of in¬ 
terest. In the postbounce phase, a neutrino leakage/heating 
scheme approximates the effects of neutrinos. We choose 
snapshots of the supernova structure beginning 3 ms after 
bounce and average the next 10 ms of evolution (sampled by 
snapshots every 1 ms) to determine the background struc¬ 
ture. This procedure smooths out most of the waves, tur¬ 
bulence, and other short-lived features present within the 
supernova core without allowing for significant evolution of 
the background structure. Our fiducial model is the A3O04 
model of A14. The rapid rotation of this model is suffi¬ 
cient to generate a rotationally dominated GW signal, but 
slow enough to be reasonably approximated by our semi- 
analytical techniques described below. 

Figure shows a density profile and propagation di¬ 
agram for the central 100 km of our fiducial model, aver¬ 
aged over the time period 3-12 ms after bounce. The inner 
~30km of the star comprise the high density PNS, which is 
surrounded by a much lower density envelope. The shock ra¬ 
dius is near r ~ 95 km during this time, and the shock does 
not strongly affect the nature of waves propagating near the 
PNS. In our models, the PNS is always stably stratified; con¬ 
vection driven by a lepton gradient does not develop until 
later times, which we do not study here. 



Figure 2. Density (p), Brunt-Vaisalii frequency (N), Lamb fre¬ 
quency {L 2 ), and angular spin frequency (f2) profiles in the inner 
100 km of the early postbounce evolution of model A3O04 of 
A14. Positive values of N'^ are denoted by the real part of the 
Brunt-Vaisala frequency (W) while negative values of corre¬ 
spond to imaginary values of the Brunt-Vaisala frequency (W). 
The PNS occupies the inner ~ 30 km of the remnant. Waves near 
the f-mode frequency have / ~ 0.8 kHz, corresponding to angular 
frequencies tu = 2iTf ~ 5 kHz, marked by the horizontal line. The 
model is made by averaging the simulation output over 10 ms, 
starting 3 ms after bounce. At these times, the shock is located 
near r ~ 95 km. 

2.2 Wave Excitation and Computation 

The sudden deceleration at core bounce excites waves which 
propagate within the PNS and surrounding material. Here, 
we semi-analytically calculate the spectrum of waves ex¬ 
cited by the bounce. We use linear and adiabatic approxi¬ 
mations for displacements from the background state, which 
we assume to be in hydrostatic equilibrium. We discuss non- 
adiabatic and non-linear effects in Section [S] We also tem¬ 
porarily ignore special/general relativistic effects and the 
impact of Goriolis and centrifugal forces, which we address 
in Section [3)^ 

Applying the approximations listed above, the lin¬ 
earized momentum equation is 

^4 = _ivjP-V<5$-5^f + f(r,t) . (1) 

aU p p 

Here, ^ is the Lagrangian displacement, p is the density, g 
is the gravitational acceleration, 5P and 5p are the Eulerian 
pressure and density perturbations, and f is the force per 
unit mass provided by the bounce. We also use the continuity 
equation 

5p+V-(p^)=0. (2) 

the adiabatic equation of state 

Sp = ^6P + , (3) 

ci g 
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and Poisson’s equation 

V^5<1> = 4nG5p . (4) 

Here, 5<£> is the gravitational potential perturbation, is 
the sound speed, N is the Brunt-Vaisala frequency, and G 
is Newton’s gravitational constant. 

Waves are excited by the force applied during the 
bounce of the inner core. For a spherically symmetric col¬ 
lapse, the strength of the force will be comparable to the 
amount of force required to halt the collapse of a shell of 
material at radius r infalling at the escape velocity. Its di¬ 
rection will be radially outwards. A rough estimate of the 
magnitude of the force per unit mass is 

f(r) ~ gr , (5) 

where g is the local gravitational acceleration in the hydro¬ 
static postbounce material, and f is the radial unit vector. 
The force peaks at the moment of the bounce, which we 
define as t = 0. The duration of the bounce is approxi¬ 
mately equal to the local dynamical time, tbounce ~ tdyn = 
\/ r^/GM{r). We approximate the time dependence of the 
bounce as a Gaussian of width equal to the dynamical time 
such that 

f(r, t) ~ f . (6) 

Note that both g and tdyn are functions of radius, so both the 
magnitude and duration of the force are strongly dependent 
on radial position. Because the gravitational acceleration g 
peaks in the core of the PNS, the forcing is concentrated in 
this region (see Figure]^. 

In a non-rotating progenitor, the collapse would be 
spherically symmetric and only radial oscillations would be 
generated. In rapidly rotating progenitors, the background 
structure is centrifugally distorted, leading to the genera¬ 
tion of axisymmetric non-radial waves. We may expect the 
degree of the non-spherical component of the force to be 
proportional to the centrifugal distortion of the collapsing 
star, e, defined such that the background density structure 
has the form 


p{r,e) = p{r)[l + e{r)Y 2 o{d)]. (7) 

Here, Y 20 is the I = 2, m — 0 spherical harmonic, and p{r) 
is the spherically averaged density profile. The centrifugal 
acceleration has magnitude ~ , where Q is the local 

spin frequency. Thus, the centrifugal distortion scales (in 
the limit e 1) as 

e ~ (H/Hdy„)" , (8) 

where Hdyn = is the local dynamical frequency. In Ap¬ 
pendix]^ we show that the perturbation in the bounce force 
per unit mass due to the centrifugal distortion is 


5f ~ A 


f^ege 

2 Y 20 f -l- r Vx V 20 

V TT 



( 9 ) 


where Vj_ is the non-radial component of the gradient, and 
A~1 parameterizes the magnitude of the force. 

With an estimate of the forcing function in hand, we can 
solve Equations |1|4| for the forced response of the fluid due 
to the bounce. It is most convenient to solve these equations 
in the frequency domain rather than the time domain. To 


do this, we decompose all perturbation variables into their 
components per unit frequency, e.g., 

^(t) = J ‘ . ( 10 ) 

Inserting this expression into Equation multiplying by 
and integrating over time, we obtain 

-a;"$„-b-V5P^-bV<5cE>,.-b<7—f = , (H) 


= 5L 


2 Y 2 QV + rV 1 Y 20 


( 12 ) 


and 

5f„ = e^tdyn . (13) 

In Equation [E each perturbed quantity is the perturba¬ 
tion per unit frequency. Similar equations can easily be de¬ 
rived from Equations |2|4| and are given in Appendix 
The Gaussian frequency dependence of the forcing term in 
Equation entails that only waves with angular frequen¬ 
cies UJ < 1/tdyn will be strongly excited at a particular lo¬ 
cation. The dynamical time is smallest at the center of the 
PNS where the density is highest, we therefore do not ex¬ 
pect significant excitation of waves with angular frequencies 
cu S> \/Gpc ~ 10 kHz. 

To solve for the frequency component response, we must 
also implement boundary conditions. At the center of the 
PNS we adopt the standard central reflective boundary con¬ 
ditions (see Appendix]^. However, in the outer regions of 
our computational domain (r ~ 250 km), there is no surface 
at which the waves will reflect [^Instead we expect the waves 
to propagate outwards until they dissipate (see Section]^. 
In the outer regions, the waves generally behave like acoustic 
(pressure) waves because their frequencies are greater than 
the local dynamical frequency. We therefore adopt a radia¬ 
tive outer boundary condition (described in Appendix 
that ensures only outgoing acoustic waves exist at the outer 
boundary. 


2.3 Wave Propagation and Characteristics 


The behavior of waves can be understood from their WKB 
dispersion relation 


k 


2 

r 




(14) 


where kr is the radial wavenumber, 

j-2 _ l{l + l)Cs 

~ ..2 


(15) 


is the Lamb frequency squared, and N'^ is the Brunt-Vaisala 
frequency squared. In regions where u> > Li and oj > N, 
waves behave like acoustic waves, while they behave like 
buoyancy (gravity) waves (not to be confused with GWs) 
where u) < Li and uu < N. Radial profiles of N and the 


^ A possible exception is the shock front, however, we expect the 
waves to be largely dissipated before they reach the shock (see 
Section]^. 
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Figure 3. Top Left: Real {Ur') and imaginary {Ui) parts of the radial wave displacement per unit frequency (Equation |A5^ for 
bounce-exited waves with / = 0.77 kHz. Top Right: Magnitude of the quadrupole moment |5Quj| (Equation |18[ l for the same waves. 
Bottom Left: Time integrated wave energy per unit radius (Equation |B5| l for waves of different frequencies. The three shown frequencies 
approximately correspond to the PNS quadrupolar (Z = 2) f-mode (/ = 0.77 kHz), the PNS I = 2 gi-mode (/ = 0.50 kHz), and an 
outgoing quasi-radial wave (/ = 1.2 kHz). Bottom Right: Strength of the excitation force f,^ (Equation |13| l. All quantities are normalized 
via their maximum values within the computational domain (see Figure]^ for dimensional values of the wave amplitude). 


quadrupolar Lamb frequency, L 2 , are shown in Figure 
Recall that typical bounce-excited waves have Z = 0 or Z = 2 
and angular frequencies near lj ~ \/Gpc ~ 5 kHz, correspond¬ 
ing to /~ 0.8 kHz. 

Near the outer boundary, L 2 and N are much smaller 
than typical wave frequencies uj, and in this limit the dis¬ 
persion relation reduces to that of acoustic waves, 


The wavelength shortens as the waves propagate outwards 
into regions with smaller sound speeds. Moreover, the WKB 


wave amplitude scales as oc and so the waves 

grow in amplitude (although their energy flux remains con¬ 
stant) as they propagate outwards. These factors cause the 
waves to become increasingly non-linear as they propagate 
outwards, i.e., they steepen into shocks. The radial group 
velocity of the waves in the outer regions is simply Vg — Cs, 
so that wave energy travels outwards at the sound speed. 

At the center of the PNS, L 2 —>■ 00 and ZV —>■ 0, so 
quadrupolar waves are evanescent for r < 5 km. In this re¬ 
gion, the response is not wave-like, but is characterized by 
the coherent fundamental mode-like oscillation of the PNS. 
The value of N'^ is large where the density gradient is large 
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at radii 5 km < r < 20 km. This region of the star can harbor 
buoyancy waves for frequencies a; < 5 kHz. 

Figure shows a plot of the radial component of the 
wave displacement per unit frequency, U, for a quadrupolar 
wave with angular frequency u ~ 5 kHz (/ ~ 0.8 kHz) as 
a function of radius. The wave contains both a real part 
(Ur) and an imaginary part (Ui), which are perpendicular 
in phase for a propagating wave and in phase for a standing 
wave. These waves reflect at the edge of the PNS (at radii of 
r ~ 20 km), so the response within the PNS is composed of 
both ingoing and outgoing waves, which interfere to produce 
a standing wave, or oscillation mode. There are no nodes 
in the radial wave function within the PNS (r < 15 km), 
therefore we refer to this mode as the fundamental mode 
(f-mode) of the PNSQ The standing waves are not totally 
reflected, and gradually leak into the surrounding material. 
For r > 50 km, the real and imaginary component of U{r) 
are perpendicular in phase, characteristic of an outwardly 
propagating acoustic wave. 

Figure also plots the strength of the forcing, f^^ (in¬ 
tegrand of Equation 131, and the time-integrated wave en¬ 


ergy per unit radius, given by Equation |B5| for waves of 
different frequencies. The forcing is localized to near the 
PNS, especially for higher frequency waves. For quadrupolar 
waves, the displacements are largest outside the inner core 
(r > 30 km), although the wave energy density is primarily 
localized to the inner ~20 km. This indicates that quadrupo¬ 
lar waves are trapped within the PNS, and only gradually 
leak out into the outer regions. For quasi-radial waves, the 
time-integrated wave energy is smaller in the PNS and larger 
in the envelope, indicating these waves are not well-trapped 
in the PNS and quickly propagate outward. We shall see in 
Sectionj^that wave energy is sharply peaked at characteris¬ 
tic frequencies that correspond to the PNS oscillation mode 
frequencies. The waves shown in Figure have frequencies 
approximately corresponding to the PNS I — 2 fundamental 
oscillation mode (f-mode, / ~ 0.8kHz), the Erst I = 2 grav¬ 
ity mode (gi-mode, / ~ 0.5kHz), and an outgoing Z = 0 
pressure wave (p-wave, / ~ 1.2 kHz)g 

Finally, Figure shows the m = 0 component of the 
wave quadrupole moment per unit frequency, 

SQ^ = j dV F5p^{r)Y2l ■ (17) 

For the Z = 2, m = 0 waves, Equation 1 1 7| reduces to 


SQu, = J dr r"‘5p„(r) . 


(18) 



Figure 4. Approximate maximum amplitudes of the radial dis¬ 
placements, ^r, associated with waves with frequencies near the 
I = 2 f-mode. We also plot the approximate non-linearity parame¬ 
ter, |fer?r|, associated with the waves. Waves are highly non-linear 
and are expected to generate shocks where |fcr?r| 1- Since the 
waves are somewhat non-linear, our linear results are only ap¬ 
proximate, and may differ from simulation results by a factor of 
order unity. Wave damping is not included in this plot, and these 
amplitudes will likely be diminished at large radii and late times. 


shown in Figure Precise wave amplitudes require an in¬ 
tegral over the response per unit frequency at a given time 
t. However, as discussed in Section the wave response 
is sharply peaked near discrete frequencies approximately 
corresponding to oscillation mode eigenfrequencies. We can 
therefore integrate the response over each of these narrow 
peaks to estimate wave amplitudes. We justify this proce¬ 
dure in Appendix |B3| The resulting amplitudes imply radial 
wave displacements of ^kilometers and velocities of a few 
percent the speed of light. These amplitudes are moderately 
non-linear, which we discuss in more detail below. 


3 WAVE DAMPING AND ROTATION 


The magnitude of the quadrupole moment is somewhat os¬ 
cillatory, but generally increases with radius. In the absence 
of wave damping, GWs are more efficiently generated as the 
waves propagate outward, however, we find below that waves 
are generally dissipated before reaching large radii. 

The approximate physical amplitude of the waves is 


^ In our models the f-mode has ui < N in much of the PNS, 
therefore it has gravity mode characteristics, which we discuss in 
more depth in Section|4.2.1| 

® We label the modes by the number of nodes in the radial dis¬ 
placement Ur within the PNS. The f-mode has no nodes, while 
the gi-mode has one node, and so on. 


3.1 Wave Damping 

The analysis presented above did not include any sources 
of wave damping. In reality the waves will damp out on 
a relatively short (~ 10 ms) timescale. Our goal here is to 
identify the primary damping mechanism and estimate the 
wave lifetime. 

Wave damping due to photon diffusion is orders of mag¬ 
nitude longer than any relevant timescales and can be ig¬ 
nored. The radiative diffusion of neutrinos, however, could 
potentially provide a significant damping mechanism. In¬ 
deed, the simulations of 012 show correlated neutrino and 
GW strain oscillations, implying that some wave energy may 
be carried away by neutrinos. However, these simulations 
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also showed very little difference between the gravitational 
waveforms with neutrino leakage turned on or off. |Ferrari| 


et al. (20031 find that neutrinos damp PNS oscillation modes 
on a neutrino diffusion time scale of ~ seconds. We therefore 
consider it unlikely that neutrinos can significantly damp 
the waves considered here on timescales as short as tens of 
milliseconds. 

GWs generated by the waves carry away wave energy 
and could potentially be an important source of wave damp¬ 
ing. As the waves propagate outward, their quadrupole mo¬ 
ment increases (see Figure]^ and so their GW energy emis¬ 
sion rate increases. However, the waves also become increas¬ 
ingly non-linear (see Figure]^. We find that waves nearly 
always become non-linear before they radiate a significant 
fraction of their energy into GWs. This is consistent with 
simulations ( |Ott||200^ |Kotake||2013p , which find that the 
energy radiated in GWs is a small fraction of the energy 
contained in fluid wave-like motions. Therefore it is unlikely 
that GW emission is a significant source of damping for most 
waves. 0 

If non-linear wave breaking occurs, the waves generate 
shocks at which point their energy is rapidly converted into 
increased entropy of the fluid where the shock forms. We 
expect the waves to non-linearly dissipate when their ampli¬ 
tude is comparable to their wavelength, i.e., when 


kr^r ~ 1 . (19) 

Figure shows an estimate of the physical displacements 
of the waves, and their degree of non-linearity. The largest 
amplitude waves (with frequencies / ~ 0.8 kHz) are moder¬ 
ately non-linear, reaching amplitudes kr^r ~ 0.5 within the 
PNS. Indeed, the simulations of A14 and K15 show evidence 
for the first harmonic of these waves in the GW spectra^ 
which is one possible outcome of non-linear wave coupling. 
In the absence of other sources of damping, the waves be¬ 
come very non-linear when r > 90 km, and will non-linearly 
break if they are able to propagate that far. We also find 
that g-modes are very non-linear within the PNS, and likely 
break and dissipate within the PNS on short timescales. 

The arguments above suggest the waves are likely dis¬ 
sipated via non-linear effects and/or turbulent damping in 
regions above the PNS (r > 40 km). An important timescale 
is the wave crossing timescale 

tcross(r) = r — , (20) 

Jo '^9 

where Vg is the wave radial group velocity, Vg = duldkr. 
Outside of the inner core (r > 30 km), Vg ~ Cs and so the 
wave crossing time for acoustic waves is approximately 



^ High frequency waves (/ ^ 2 kHz) may radiate much of their 
energy in GWs due to the dependence of the GW energy 
emission rate. However, since these waves are weakly excited at 
bounce, they contain little energy and cannot generate a strong 
GW signature. 

^ These peaks are more prominent in postbounce spectra, i.e., 
spectra where the bounce is windowed out. This allows the non¬ 
linear harmonic peak to be separated from the broad spectrum 
of GWs contributing to the bounce signal. 



/ (kHz) 


Figure 5. Wave crossing timescale tcross and wave leakage 
timescale tieak across the inner 50 km of the postbounce super¬ 
nova structure as a function of wave frequency for 1 = 2 and Z = 0 
waves. Each peak in tieak corresponds to a PNS oscillation mode; 
the peak at / ~ 0.77 kHz is the 1 = 2 f-mode. Quasi-radial (Z = 0) 
waves have tcross ~ tieak because these waves propagate outward 
without significant reflection. 


The last equality arises because the regions of wave propaga¬ 
tion are approximately virialized, so that the sound crossing 
time is comparable to the dynamical time. Low frequency 
waves are buoyancy waves in the PNS, and these waves have 
much larger wave crossing times because the radial group ve¬ 
locity of buoyancy waves is much smaller than that of sound 
waves. Figure shows tcross evaluated at r = 50 km as a 
function of wave frequency. Wave crossing times over this 
region range are ~ 1 ms for high frequency waves. 

However, tcross (r) is not always a good estimate for a 
wave damping timescale because waves can reflect off the 
PNS surface and become trapped within the PNS. Low fre¬ 
quency waves only gradually tunnel through the overlying 
evanescent region (cf. Figure]^, and therefore a more rel¬ 
evant timescale is the wave leakage timescale tieak (r). This 
timescale reflects the rate at which wave energy leaks out 
of regions below r, and is calculated via Equation |B13| in 
Appendix |B2| 

Figure 1^ plots tieak, evaluated at r = 50 km, as 
a function of wave frequency. For high frequency waves 
(/ > IkHz), tcross — tieak because these waves are not 
trapped within the PNS. They should therefore damp out 
on timescales of milliseconds. Lower frequency waves exhibit 
a series of peaks in the value of tieak. These peaks corre¬ 
spond approximately to PNS oscillation modes, for which 
the waves are mostly reflected at the PNS surface. Waves 
with frequencies / ~ 0.77 kHz, corresponding to the PNS 
f-mode, leak outwards on timescales of tieak ~ 10 ms. This 
timescale is quite similar to the damping times found in sim¬ 
ulations (e.g., 012, A14). We therefore conclude that the 
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lifetime of these waves is well approximated by tieak- Waves 
at this frequency are likely damped via turbulent/non-linear 
effects only after they are able to leak into the outer enve¬ 
lope (r > 30 km). Waves at lower frequencies corresponding 
to gravity modes within the PNS may have much longer life 
times, although we shall see in Section that their contri¬ 
bution to the GW spectrum is quite small. 


3.2 Rotation 


The most obvious impact of rotation is to provide centrifugal 
support for the PNS and surrounding material. As the an¬ 
gular spin frequency Q increases, centrifugal support causes 
the postbounce structure to be less compact (i.e., lower cen¬ 
tral densities and larger PNS radii) and more oblate. The 
change in background structure is well-captured by the sim¬ 
ulations A14 used to generate our background structures. 
We have also attempted to approximately account for the 
effect of the oblateness in the strength of the forcing that 
excites the waves (see Section]^. Because the centrifugally 
supported PNSs are less compact, their corresponding dy¬ 
namical frequencies fldyn and mode frequencies are gen¬ 
erally lower. Including only the effect of rotation on the 
background structure, we would expect the frequencies of all 
axisymmetric modes to decrease with increasing rotation. 

However, in the wave calculations presented in Ap¬ 
pendix 1^ we have ignored the Coriolis and centrifugal 
terms in the momentum equations. These terms become 
very important when the oscillation mode frequency / be¬ 
comes comparable to 2/spin. We are primarily interested 
in waves with frequencies / < 1 kHz, whereas the inner 
core spin frequencies of the background models are of order 
2/apin ~ 0.6 kHz. Rotation therefore has a strong effect on 
the wave dynamics and may significantly change our results. 
Details on the influence of rapid rotation on the oscillation 


modes of NSs can be found in Bildsten et al. (19961, Dim- 
[melmeier et al.| ( |2006[ ), and [Passamonti et al.| ( [2009[ ). 

We can attempt to predict the influence of rotation 
based on perturbation theory. For the axisymmetric oscil¬ 
lation modes of interest, the first order change in mode fre¬ 
quency (proportional to H) vanishes. Second order correc¬ 
tions include the Coriolis and centrifugal forces, the non- 
spherical background structure, and spin-induced coupling 
between oscillation modes. For axisymmetric I = 2 f-modes. 


rotation has only a small effect on the mode frequency (Dim- 
[melmeier et al.||2006[ ). For low frequency axisymmetric g- 
modes, however, rotation typically increases axisymmetric 
mode frequencies (see Bildsten et al.jl9'^ Lee & Saio 19971. 
The f-mode with / ~ 0.8 kHz in Figureis somewhat mixed 
in character. It behaves like an f-mode at r <10 km where 
it is evanescent, but behaves like a g-mode in the range 
10 km < r < 20 km. Its mixed character makes is difficult to 
easily predict the effect of rotation on its frequency. 

An additional complication is that rotation couples 
spherical harmonics of Yi^rn and Yi± 2 ,m- Axisymmetric 1=2, 
m = 0 waves of interest couple to both 1 = 0 (radial) waves 
and 1=4 waves. Consequently both 1=0 and 1=4 waves will 
obtain quadrupole components that allow them to generate 
GWs. Rotation also induces mode mixing between modes of 
the same degree I, e.g., between the 1=2 PNS f-mode and g- 
modes. Rotational mixing between the modes prevents their 


frequencies from crossing (see Section 4 of Fuller et al. 20141, 


and the modes instead undergo an “avoided” crossing in 
which they exchange character. During the avoided cross¬ 
ing, the mode frequency separation is approximately con¬ 
stant, and the modes are superpositions of the unperturbed 
modes, giving them a hybrid mode character. We revisit the 
observational consequence of rotational mode coupling in 
Section 14.2.11 

Finally, rotation introduces inertial waves/modes, 
which are restored by the Coriolis force. In uniformly ro¬ 
tating bodies, inertial modes generally exist within the an¬ 
gular frequency range uj < 2H. In our most rapidly rotating 
models, the f-modes and g-modes lie in the range uj < 2Hc 
(where Qc is the peak angular velocity in the postbounce 
model). Therefore inertial modes could potentially influence 
the wave spectrum, either directly (by producing a peak in 
the GW spectrum) or indirectly (by affecting the frequency 
of the f-mode). 


3.3 Special and General Relativity 


In this work, we largely ignore the effects of special rel¬ 
ativity (SR) and general relativity (GR). Although rele¬ 
vant for the PNS, they greatly complicate the analysis. 
The effects of GR on NS oscillation modes have been stud- 


Cutler & LindblomI 

19921 Andersson 19981 Lockitch et al. 

2001[ 20031 Boutlou 

cos & Nollert 20071 Gaertig & Kokko- 


can anticipate GR to affect mode frequencies by at most 
0\2GM{ r)/{rF)] < 10% at all radii within our model. 
However, Dimmelmeier et al.|(2002l compared Newtonian 


and conformally flat GR CG simulations, finding significant 
differences in postbounce structure and GW spectra. Since 
our background structures are generated from GR simula¬ 
tions, we expect that GR effects on wave dynamics will be 
smaller than the effects of rapid rotation. 

SR effects also become important when fluid motions 
approach the speed of light. The velocity associated with a 
perturbation with = 2 km at a frequency / = 1 kHz is 
Sv « 1.2 X 10® cm/s « 4 X 10“®c. Therefore the wave ampli¬ 
tudes predicted by our calculations (see Figure]^ generate 
fluid velocities far below speed of light, and so SR corrections 
are small. 


4 GRAVITATIONAL WAVE SIGNATURES 


GW and and neutrino emission are the only way of directly 
observing the core dynamics of CG SNe (Ott 20091. Here 
we attempt to quantify the GW signatures produced by the 
bounce-excited oscillations and compare them with simula¬ 
tion results. 


4.1 Gravitational Wave Spectrum 

We now turn our attention to GW wave emission induced 
by the fluid waves. The time-integrated GW energy emitted 
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per unit frequency is 

27r J dt Egw,ll 


_ dEaw 
~ df 
2G 6ixn 

= \SQ^ 


( 22 ) 


where SQui is the quadrupole moment per unit frequency 


(calculated from Equation 171 and the second line is from 


012. The GW energy corresponds to a characteristic dimen¬ 
sionless wave strain ([Flanagan fc Hughes||1998| 


hchar — 


2 G dEGw 
TT^ (PD^ df 


(23) 


We use the fiducial distance D — lOkpc in our presented 
results. 

Caution must be used when evaluating Equation |22| 
Although the wave frequency is constant in radius, the 
quadrupole moment \SQui\ generally increases at larger radii 
(see Figure]^. The GW energy flux is therefore dependent 
on which radius we choose to evaluate Moreover, the 

total energy emitted by the GWs could be larger than the 
wave energy, especially for high frequency waves. This un- 
physicality reflects the fact that we have not taken GW emis¬ 
sion into account in the fluid oscillation equations; in reality 
the waves are attenuated as they emit GWs. 

In what follows, we calculate GW energies and ampli¬ 
tudes with SQui evaluated at row = 30 km. This radius is a 
good choice as long as the fluid waves damp out at radii just 
above row - If they are able to propagate to larger radii, they 
will obtain larger quadrupole moments (see Figure and 
may emit more energy in GWs. Therefore, the energy fluxes 
we calculate should be viewed as order of magnitude es¬ 
timates, and only full non-linear hydrodynamic simulations 
can yield quantitatively reliable predictions. The frequencies 
of the peaks in the GW spectrum are not strongly affected 
by our choice of raw because these peaks are primarily de¬ 
termined by the values of the PNS mode frequencies. 

Figure shows a plot of the total energy radiated in 
GWs per unit frequency and the associated characteristic 
wave strains. For frequencies / < 1 kHz, the GW energy 
is sharply peaked around characteristic frequencies. These 
frequencies essentially correspond to the oscillation mode 
frequencies of the PNS. The widths of the peaks are the 
inverse of the mode lifetimes, i.e., the timescale on which 
the mode energy is able to leak out of the PNS, tieak (see 
Equation B131. There are no sharp peaks at high frequencies 


/ > 1 kHz because these frequencies correspond to pressure 
waves (which are not reflected at the edge of the PNS) that 
quickly propagate outwards on a wave crossing time tcross. 

The peak at / ~ 0.8 kHz corresponds to the 1=2 axisym- 
metric PNS oscillation mode (see Section [2.3| ). This mode 
contains more energy than any other because it couples well 
(in both physical and frequency space) to the forcing pro¬ 
duced at bounce. We therefore confirm the hypothesis of 012 
that the peak centered at / ~ 0.8 kHz in their GW spectra 
is generated by the axisymmetric quadrupolar PNS f-mode. 
The f-mode is expected to dominate the early postbounce 
GW signature of rapidly rotating CG SNe. The peaks at 
lower frequencies correspond to 1=2 PNS g-modes; the hrst 
is the gi-mode at / ~ 0.5 kHz. We reiterate that rotational 
effects are important and may substantially alter mode fre- 



Figure 6. Top: GW energy spectrum, dFow/df, due to bounce- 
excited oscillations, calculated using our A3O04 model (from 
A14). We have plotted contributions from 1 = 2, 1 = 4, and Z = 0 
waves. The 1 = 0 and I = 4 waves emit GWs because they gain a 
quadrupole moment due to the aspherical background structure, 
and their energy spectrum is only an order of magnitude estimate 
(see text). The broad spectrum of GWs at higher frequencies is 
created by outgoing pressure waves that create the bounce signal. 
This signal is likely to be overestimated, as these waves quickly 
dissipate via shock formation before they can radiate GWs. Bot¬ 
tom: Dimensionless characteristic wave strain (Equation |23| l for 
the same model, computed at D = 10 kpc. We have labeled peaks 
corresponding to some of the PNS oscillation modes. 


quencies, especially in the low frequency part of the spec¬ 
trum. 

We have also attempted to calculate the GW spectra 
due to quasi-radial {I = 0) and I = 4 waves, which obtain 
quadrupole moments due to rotational mixing with 1 = 2 
modes (see Section 3.21 and due to the centrifugally dis¬ 
torted background structure. To estimate GW spectra for 
1 = 0 waves, we remove the e = (fl/fldyn)^ dependence of 
the forcing term (see Equation |13| ) because centrifugal dis¬ 
tortion is not needed to excite radial waves. To calculate 
their quadrupole moment, we multiply the right hand side 
of Equation by e, which approximately accounts for the 
quadrupole moment of the background structure and allows 
the 1=0 waves t o em it GW. For Z=4 waves, we replace e with 
to estimate the reduced strength of the 


in Equation 


13 


wave excitation, and use the same procedure as Z = 0 waves 
to calculate an approximate quadrupole moment. This pro¬ 
cedure is rudimentary and should not be expected to yield 
accurate quantitative predictions for the energy radiated by 
I = 0 and I = 4 waves, although it can be used for a quali¬ 
tative understanding. 
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The GW spectrum produced by the 1 = 0 waves is a 
smooth continuum rather than being peaked at mode fre¬ 
quencies. The reason is that 1=0 waves are not well reflected 
from the PNS edge, and so energy in 1=0 waves leaks out of 
the PNS on a wave crossing (dynamical) timescale. More¬ 
over, 1=0 g-modes do not exist, instead low frequency 1=0 
waves are evanescent in the PNS when co < N. This is in 
stark contrast to the 1 = 2 waves, which can be trapped in 
regions with cu < to form oscillation modes. Instead, the 
force exherted by the bounce is transferred to 1=0 waves of 
a broad range in frequencies, which quickly travel outward 
and steepen into shocks; it is this process which generates 
the outgoing shock created by bounce. Although the 1 = 0 
waves are important for the GW spectrum of Figure for 
/ > 1 kHz, their GW strain peaks near bounce and con¬ 
tributes primarily to the bounce signal (see Figure and 
Section]^. The same is true for higher frequency (/ > 1 kHz) 
quadrupolar waves. After bounce, the high frequency waves 
are quickly dissipated by shocks, and the 1=2 modes dom¬ 
inate the GW spectrum. This idea is consistent with the 
results of K15, who find the GW spectrum is more strongly 
peaked around mode frequencies when the bounce is win¬ 
dowed out. 

The 1=4 waves are less efficiently excited by the bounce 
than 1=2 waves, but the response is similarly peaked around 
mode frequencies. The largest peak is centered around the 
1=4 f-mode at / ~ 1 kHz, although we expect this mode to 
radiate considerably less GW energy than the 1=2 f-mode. 
Nonetheless, given our rudimentary methods, we speculate 
that the 1=4 f-mode may be detectable, especially for very 
rapidly rotating progenitors. 


4.2 Comparison with Non-linear Simulations 

We now compare our semi-analytical results with the sim¬ 
ulation results of 012 and A14. The GW energy spectra of 
these simulations generally contain a few distinguishing fea¬ 
tures: 

1. A prominent peak of maximum GW energy in the range 
0.7kHz < /max < 0.8 kHz. 

2. A broad spectrum of GW energy at frequencies / < 2 kHz. 

The main peak at /max ~ 0.75 kHz is due to the 1=2, 
m = 0 f-mode of the PNS, as speculated by 012. However, 
both the frequency and wave function of this mode are likely 
influenced by rotational interaction with other modes (see 
Section [3. 2| and discussion below). The gi-mode may be re¬ 
sponsible for peaks near / ~ 0.5 kHz, and the 1=4 f-mode 
may produce a peak at / ~ 0.95 kHz. Our linear calcula¬ 
tions do not easily account for any peaks at / > 1 kHz. We 
speculate that peaks near / ~ 1.5 kHz are due to the first 
harmonic of /max, generated due to the non-linearity of the 
f-mode responsible for /max (see Appendix |B3[ ). 

The broad spectrum of GW energy (visible as broad 
peaks near / ~ 1.2 kHz in Figure and the top panel of 
Figure]^ is produced by both quasi-radial and quadrupolar 
waves which are not efficiently reflected and quickly propa¬ 
gate out of the PNS. However, we caution that our methods 
may overpredict the GW signal from these waves, as they 
quickly steepen into shocks before generating GWs. This 
may account for the lack of a very broad peak at / ~ 1.2 kHz 
in the rapidly rotating simulations of 012 and A14 (lower 
panel of Figure [^. 


Our calculations predict total GW energy outputs and 
wave strains roughly consistent with the results of 012 and 
A14 when we use a forcing strength of A ~ 2 (Equation |13[ ), 
although there are significant uncertainties in calculating 
the GW spectrum. Here, we claim only that our method 
produces a sensible order of magnitude estimate for GW 
energies, and that it provides a physical explanation for some 
of the features in the GW spectra from rotating core collapse 
simulations. 

Finally, we comment on the widths of the GW spec¬ 
tral peaks. As discussed in Section]^ the damping timescale 
for the oscillation modes is determined by the wave leak¬ 
age timescale into the envelope. For the 1 = 2 f-mode, this 
leakage time is ~ 10 ms (in good agreement with the GW 
decay timescale seen in 012 and A14), corresponding to the 
width of A/ ~ 0.1 kHz for the peak at /max. The g-mode 
peaks in Figure]^ are narrower on account of the long leak¬ 
age timescale for the g-modes, but their widths are under¬ 
estimated since the g-modes may be damped via non-linear 
processes or modified by the background structural evolu¬ 
tion. 


4-2.i Rotation 


A14 computed GW energy spectra of rapidly rotating pro¬ 
genitors as a function of the rotation rate. To compare with 
their results, we generate background structures and pre¬ 
dicted GW spectra from the simulations of A14 for different 
rotation rates, using the same techniques described in the 
preceding sections. 

Figure [^displays our computed GW spectra (top panel) 
for the rapidly rotating models of A14, in addition to the 
GW spectra obtained from the simulations themselves (bot¬ 
tom panel). As the rotation rate increases, the PNS and 
surrounding material have more centrifugal support and be¬ 
come less compact, decreasing their dynamical frequency. 
Lamb frequencies L 2 , and Brunt-Vaisala frequencies N. For 
this reason, the prominent peaks computed from our semi- 
analytical analysis shift to lower frequencies with increas¬ 
ing rotation rate. They also decrease in power due to the 
strong dependence of Eqw on frequency (see Equation 221. 
Rapidly rotating models also contain significant power in a 
broad peak centered near / ~ 1.2 kHz. This power is cre¬ 
ated by pressure waves which quickly leak into the envelope, 
although this power may be overestimated (see discussion 
above). 

A monotonic trend of decreasing /max with spin fre¬ 
quency is not observed in A14. Instead, their results indi¬ 
cate a more complex dependence in which /max increases in 
frequency at moderate rotation rates and decreases at very 
high rotation rates, although over most of the parameter 
space 0.7 kHz < /max < 0.8 kHz. 

We attribute the difference to our neglect of the ef¬ 
fects of the Goriolis and centrifugal forces on the wave dy¬ 
namics, although GR effects could have some influence as 
well. As discussed in Section [3.2| the f-mode is really a hy¬ 
brid f-mode/g-mode that propagates near the sharp density 
gradient at the interface between the PNS and surround¬ 
ing material. Due to its g-mode component, the Coriolis 


force will tend to increase its frequency (Bildsten et al. 
|1996p in rapidly rotating models. We speculate that the in¬ 
creasing impact of the Coriolis force with rotation largely 
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Figure 7. Top: Our computed GW spectra for the sequence of 
models A3002-A3008 of A14. Larger model numbers are more 
rapidly rotating. The spectrum of each model has been scaled to 
its maximum amplitude and shifted vertically for clarity. The f- 
mode corresponds to the sharp peak in each curve, and is marked 
with a circle symbol. The broad humps of power at higher fre¬ 
quencies are due to outgoing pressure waves which create the 
bounce signal, although their contribution to the GW spectrum 
is likely overestimated (see text). Bottom: GW energy spectra 
from the same simulations (A14) used to generate the background 
models for our calculations. The peak which we speculate to be 
generated by the quadrupolar f-mode has been marked with a cir¬ 
cle symbol, although its location is ambiguous in slowly rotating 
models. The mismatch in frequency of the f-mode peak between 
top and bottom panels is most likely due to the strong effect of 
rotation on the wave dynamics (see text), which is not included 
our calculations. 

cancels the effect of the less compact structure, causing 
0.7 kHz < /f-mode < 0.8 kHz over a large range of spin fre¬ 
quencies. A fairly weak dependence on spin frequency is also 
seen for axisymmetric Z = 2 f-modes in IDimmelmeier et al.l 
( 2 ^. 

Rotational mixing between the Z=2 f-mode and gi-mode 
may also affect the frequency behavior of the GW spectra. 
This mixing can generate mode “repulsion” between the f- 
mode and g-mode during an avoided crossing, which may 
increase the value of /f-mode at some spin rates. Moreover, 
in very rapidly rotating models for which Wf-mode 2f2c, 
mixing with inertial modes may also influence the wave spec¬ 
trum. 

We also note that 012 find an oscillation in the cen¬ 
tral density pc of their simulations with frequency /max. If 
the oscillation at /max is generated by a purely quadrupolar 
mode, we should expect a negligible oscillation in pc because 


all perturbations vanish at the center for modes with I ^ 2. 
However, as described above, rotation induces mixing be¬ 
tween waves of I and Z ± 2. In rapidly rotating models the 
quadrupole waves thus obtain a radial component which will 
generate oscillations in pc. 

Properly incorporating the influence of rotation is be¬ 
yond the scope of this paper. We hope that future studies 
with more rigorous implementations of rotational effects will 
improve our understanding of their influence on the GW 
spectrum. Unfortunately, the rather complex dependence of 
/max on rotation frequency may complicate interpretations 
of a GW signal from a nearby GC SNe. Simulation results 
(which naturally include effects of rotation, non-linearity, 
etc.) such as those of A14 (and the upcoming study by K15) 
may therefore be vital for interpretation of observational re¬ 
sults. 

4.3 Neutrino Signatures 

Neutrinos offer another opportunity to observe the interior 
dynamics of CC SNe. Our calculations do not include any ef¬ 
fects of neutrinos, but we may speculate on the neutrino sig¬ 
nature of bounce-excited oscillations. The oscillations gen¬ 
erate perturbations in density in temperature which may 
generate oscillations in neutrino generation rates. The oscil¬ 
lations may also perturb the effective location of the neu- 
trinosphere. Both these effects may generate oscillations in 
the observed neutrino flux, with identical frequencies to the 
peaks in the GR spectrum. Indeed, 012 find oscillations in 
the neutrino flux that are correlated with the GR signal. We 
do not investigate the physics of these oscillations in detail, 
but merely speculate that we expect a peak in the Fourier 
transform of a neutrino light curve at 0.7 kHz < / < 0.8 kHz. 


5 DISCUSSION AND CONCLUSIONS 

We have examined the physics of wave generation, propaga¬ 
tion, and dissipation in rapidly rotating core-collapse super¬ 
novae (CC SNe). Our linear, semi-analytic methods comple¬ 
ment recent hydrodynamic simulations of rotating CC SNe 
(012 and A14), and they shed light on the basic physics 
at play. The quadrupolar distortion of the centrifugally dis¬ 
torted PNS inner core generates axisymmetric quadrupolar 
(Z = 2) waves at bounce. Because of their quadrupolar na¬ 
ture, these waves can generate a strong gravitational wave 
(GW) signal that could be detected by advanced LIGO for 
a galactic CC SN (012, A14). 

The GW signal is largely determined by two factors: 
the strength of the wave excitation provided by core bounce, 
and the ensuing propagation/dissipation of the waves. The 
strength of the force (per unit mass) provided by the bounce 
is roughly proportional to the local gravitational accelera¬ 
tion (see Equation]^, therefore, the wave excitation is spa¬ 
tially confined to the inner core of the PNS (r < 30 km). The 
duration of the force is approximately the dynamical time 
of the inner core of the PNS, tdyn ~ 1 ms. Hence, wave ex¬ 
citation is peaked around frequencies / ~ ~ 1 kHz. For 

quadrupolar waves, the strength of the forcing is also pro¬ 
portional to the centrifugal distortion, e, of the inner core 
of the PNS, and therefore more rapidly rotating progenitors 
typically transfer more power into quadrupolar oscillations. 
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The behavior of quasi-radial (I — 0) and quadrupolar 
(Z = 2) bounce excited waves is fundamentally different. Un¬ 
like quadrupolar waves, the quasi-radial waves are not re¬ 
flected at the edge of the PNS and quickly propagate into the 
overlying envelope where they generate the bounce-induced 
shock. Although the GW bounce signal (see Figure]^ con¬ 
tains a contribution from the quasi-radial bounce of the PNS 
(which has a quadrupolar component because of the cen¬ 
trifugal distortion of the PNS), quasi-radial waves contribute 
little to the GW ring down signal after bounce. 

In contrast, quadrupolar waves propagate as buoyancy 
waves within the PNS in regions where {lj < N). These 
waves are efficiently reflected at the edge of the PNS (r ~ 
30 km), and become trapped within the PNS. In-going and 
outgoing waves interfere to create standing waves, or oscil¬ 
lation “modes”, of the PNS. Gonsequently, the quadrupo¬ 
lar wave energy is peaked around characteristic “mode” 
frequencies. The wave energy leaks out of the PNS on 
timescales of ~ 10 ms, propagating outwards into the en¬ 
velope as acoustic waves, and eventually dissipating. 

The largest early postbounce {t — tbounce < 20 ms) GW 
signal is produced by the fundamental quadrupolar oscil¬ 
lation “mode” of the PNS, whose frequency is typically 
/ ~ 0.75 kHz. This f-mode is essentially a surface wave of 
the PNS, and owes its existence to the sharp density gradi¬ 
ent at the edge of the PNS. However, we emphasize that the 
f-mode has some buoyancy wave (g-mode) character because 
it propagates in a region where oj < N. Moreover, the f-mode 
energy leaks into the low density surrounding envelope, and 
its properties are not accurately captured by a calculation 
of oscillation modes in an isolated PNS or NS. Weaker post¬ 
bounce GW signals may be produced by higher order I = 2 
g-modes or the I = 4 f-mode (like the quasi-radial I = 0 waves, 
the I = 4 waves obtain a quadrupole moment due to the cen¬ 
trifugal distortion of the background structure). Finally, the 
first harmonic of the f-mode, at a frequency / ~ 1.5 kHz, 
may appear in GW spectra due to the somewhat non-linear 
amplitude of the f-mode. 

The greatest uncertainty and source of error in our anal¬ 
ysis is the neglect of the effects of Goriolis and centrifugal 
forces on the wave dynamics. Indeed, these forces are quite 
important since typical angular spin frequencies are compa¬ 
rable to the wave frequencies. Our calculations (incorrectly) 
predict that the quadrupolar f-mode frequency monoton- 
ically decreases with increasing progenitor spin frequency 
due to the less compact background structure. In contrast, 
simulations (A14, K15) show only a weak dependence of the 
f-mode frequency on the spin frequency. This suggests that 
rotation helps increase the f-mode frequency and largely off¬ 
sets the frequency decrease due to lower average density in 
more rapidly rotating models. We speculate that the Coriolis 
force increases the f-mode frequency, as it does to g-modes 


in rapidly rotating stars (Bildsten et al. 1996 Lee & Saio 


1997). In any case, the bounce signal likely provides a bet¬ 


ter indication of progenitor spin frequency than the f-mode 
frequency (A14). 

Nonetheless, the detection of both the bounce signal 
and the f-mode would provide useful constraints on the 
structure and spin frequency of the inner core. The fre¬ 
quency of the f-mode varies little with rotation rate and 
depends primarily upon the equation of state (A14, K15), 
while the bounce signal is sensitive to both equation of state 


and spin frequency. Measuring /f_mode could thus break a 
partial degeneracy between the compactness of the inner 
core and its spin frequency. Additionally, for slow-moderate 
rotation rates, the amplitude |af_mode| of the quadrupolar 
f-mode signal is roughly proportional to the square of the 
core spin frequency, A measurement of |af_inode|, cou¬ 
pled with an accurate distance measurement (likely to be 
fairly well-determined for a galactic GG SN), will thus pro¬ 
vide an additional constraint on inner core spin frequency. 

The simple analysis presented here has shed some light 
on wave dynamics within rapidly rotating CC SNe, and how 
the basic properties of the progenitor contribute to the GW 
spectrum produced by waves excited at bounce. However, 
our analysis is not sufficient to accurately calculate GW 
spectra suitable for comparison with observations. We en¬ 
courage further simulations of rapidly rotating CG SNe that 
generate GW spectra over reasonable ranges in parameter 
space (e.g., using a range of rotation frequencies, differen¬ 
tial rotation profiles, equations of state, neutrino approxima¬ 
tions, etc.) complementary to those already performed (e.g., 
IDimmelmeier et al.| 120081 [Abdikamalov et al.||2010l 012, 
A14). These simulations will provide templates for compari¬ 
son with the GWs observed by advanced LIGO in the event 
of a galactic CC SN, and they will therefore be an important 
tool for understanding the extreme physics of these explo¬ 
sions. 
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APPENDIX A: OSCILLATION EQUATIONS 
AND BOUNDARY CONDITIONS 


A1 Force Exerted by Bounce 

We begin by estimating the magnitude of the force exerted 
on the fluid during core bounce. The force on a spherical 
shell is approximately equal to the force required to halt the 
free fall of the shell. Therefore we require 


Then, using the form of f from equation we find 

f = g . 


(Al) 


(A2) 


In a centrifugally distorted star, both the magnitude 
and direction of the force in equation A2 are perturbed (al¬ 
though we will ignore perturbations in the time dependence 
because its functional form is only approximate), such that 
the centrifugal perturbation to the force is 


/2 ^ 

Sg ^ r- 
— r J- (5n 

^ TV 

. 9 


(A3) 


Here, Sg is the perturbation in gravitational acceleration 
due to the centrifugal distortion, and Sn is the perturba¬ 
tion in surface normal to each centrifugally distorted shell. 
To first order, the centrifugal distortion perturbs the lo¬ 
cation of each spherical shell, located at radius r, by an 
amount 5r = —trY 2 o, with e ~ (B/Bdyn)^. The pertur¬ 
bation in the gravitational acceleration is approximately 
5g/g ~ —25rjr, while the perturbation in the surface nor¬ 
mal is Jn = — V_L(5r. Then the perturbed bounce force due 
to centrifugal distortion is 


e^e l‘Adyn) 


2Y2of+ rVxY2( 


(A4) 


and A is a constant of order unity that parameterizes the 
strength of the force. 


A2 Oscillation Equations 

Next we describe our method of solving the linearized forced 
oscillation Equations [T]|^ As described in the text, we de¬ 
compose the perturbed fluid variables into their frequency 
components. For brevity, we drop the lu subscript used in 
the text, but it should be understood that all perturbed 
variables are the fluid response per unit frequency. We have 
decomposed the force into spherical harmonics, and do the 
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same for the fluid response, as each spherical harmonic com¬ 
ponent will couple only to the associated forcing term. The 
radial component of the Lagrangian displacement is 

Ur) = U{r)Yi^{e,cj,)r (A5) 

while the horizontal component is 

iF^) = V{r)rVYi^{9,<t>) . (A6) 


In a rotating star there will also exist a toroidal component 
to the horizontal displacement, however computation of this 
term requires the inclusion of Coriolis and centrifugal forces 
and complicates the procedure. We proceed without includ¬ 
ing rotational effects, with the understanding that waves, 
especially at low frequencies, will be significantly altered by 
the effects of rotation. 

We define the enthalpy perturbation 


^(r) = -b h-l>(r) - r hf Y 20 


(A7) 


where 


5f = 




■egfdyne . 


(A8) 


After integrating over angle, the horizontal components 
of the momentum equation (Equation 111 yield V{r) = 

Upon angular integration, the radial component of 
Equation becomes 

^5P+ -^5P + p5g + p{N^ -uj^)U = 2pSf . (A9) 

OT Cc 


We have dropped the (r) dependence of the variables for con¬ 
venience. Transforming to the frequency domain, the conti¬ 
nuity equation becomes 


9 rr 


r c: 


U + 

5f . 


1(1 + 1) 


SP l{l + l) 


5$ 


(AlO) 


Finally, Poisson’s equation can be written as the two first 
order differential equations 




(All) 


_9 

dr 


5g -b -5g 
r 


1(1 + 1 ) 


— AnGp 


-^5P + 
pci 



= 0 . 

(A12) 


It is important to understand that Equations |A9|A12| are 
complex, and each perturbed variable has both a real and 
imaginary component. The forcing term 5f is purely real 
and generates differences between the real and imaginary 
components responsible for energy transport. 

To solve Equations |A9|A12| we need eight boundary 
conditions as required for the eight variables composed of 
the real and imaginary parts of [7, 5P, <5$, and 5g. The four 
inner boundary conditions are the usual relations 


and 


(A13) 

(jj^r 

5g = U<1> . (A14) 


As justified in the text, we require a radiative outer 
boundary condition. To do this, we use the WKB approx¬ 
imation for the variables at the outer boundary, such that 
U oc with kr — ±uj/Cs. Siuce we have defined the 

time dependence of the waves to be proportional to 


the positive value of kr corresponds to an outgoing wave. In 


the WKB limit. Equation |A10| is 


ikrll -1- T^SP — 0 , 

pci 

(A15) 

or 


U = ^ 5P . 

UjpCs 

(A16) 

Similarly, Equation |A12| is approximately 


UCs 

(A17) 


Equations I A16| and I A17| constitute our four outer boundary 
conditions. 


APPENDIX B: WAVE DYNAMICS 
B1 Work and Energy Flux 


As a check on our numerical calculations and our interpreta¬ 
tion of their results, it is helpful to understand the wave en¬ 
ergetics. The basic idea is that the forcing term in Equation 
|A10| imparts energy into the waves, and this energy is even¬ 
tually carried away by the outwardly propagating waves. 
The rate of energy change in the volume interior to radius 
r due to the waves (in the absence of any forcing) is 


E = -r^ J dA{5P + p5^)Svr - j dVS^^Sp . (Bl) 

where dA is a surface area element of the spherical surface at 
radius r. The first term on the right hand side is the energy 
flux through the surface and the second term is the gravita¬ 
tional work done by the volume. To evaluate this expression, 
we decompose the perturbations into the frequency domain 
as in Equation 1 10[ We then integrate over angle and time to 
obtain the net energy outflow per unit frequency. The result 
is 

J dtEui = 2ujr^ — Ur,ui'^i,Lj) 

+ 2ujr^ [ dr(5pr^tj5^i^tj — 5pi^^5^r,ui) ■ (B2) 


Here, the r and i subscripts denote the real and imaginary 
subscripts of each variable, respectively, and we have rein¬ 
troduced the uj subscript to denote the response per unit 
frequency. The total energy flux (including waves of all fre¬ 
quencies) is found by integrating Equation B2 over all fre¬ 
quencies. 

We can compare the energy flux of Equation |B2| with 
the energy deposited into the waves by the forcing term. The 
work done per unit time by the force is 


E = 


dVpf ■ hv . 


(B3) 


Once again decomposing the response into its spherical har¬ 
monic and frequency components, and integrating over all 
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time, we find the total work done by the force per unit fre¬ 
quency, 

J dtEui = 2a; drpr^fu, ■ „ . (B4) 

Since there are no other wave driving or damping terms, the 
total work done by the force (Equation |B4[ | should be iden¬ 
tical to the total energy outflow of Equation |B2| at all radii 
r and at all angular frequencies a;. To check our numerics 
we verify that this is indeed the case. 

Finally, we can compute the time integrated energy per 
unit radius per unit frequency. The physical meaning of this 
quantity is the energy density of a wave packet of angular 
frequency uj, weighted by the amount of time it spends at 
any location. The result is 

/ =2uj^pr^[\U^f+l{l + l)\V:.f] , (B5) 

which is plotted in Figure The waves spend most of their 
time in the core, so their time-weighted energy density is 
largest there. They travel through the outer regions rela¬ 
tively rapidly, without significant reflection, and their en¬ 
ergy density is small there. In this case their energy density 
should be related to the total energy outflow rate via 

J dtE^ J dt cs^ , (B6) 

because Cs is the group velocity of the waves through this re¬ 
gion. We verify that Equation |B6| is approximately satisfied 
in the outer layers of our computational domain. 


B2 Timescales 

It is important to understand the timescales involved for 
bounce-excited waves in PNSs. The waves of interest have 
oscillation periods comparable to both the spin period and 
dynamical timescale of the PNS: 

P ~ Ps ~ Pdyn ~ 1 ms . (B7) 

In contrast, the PNS and surrounding envelope evolve over 
longer timescales set by the mass accretion rate, 

Pevol ^ 50 ms . (B8) 

The waves have plenty of time to oscillate before they are al¬ 
tered by the evolution of the supernova. However, we should 
be skeptical of any wave timescales longer than Pevoi, as 
any processes acting on these timescales will likely be irrele¬ 
vant compared to the dynamical evolution of the background 
structure. 

To understand timescales associated with the waves, it 
is helpful to construct a toy problem. We consider the same 
supernova background structure shown in Figure Rather 
than calculate a forced wave solution as we do in Section 
we consider the properties of steady oscillations ocurring at 
an angular frequency a;. To do this, we solve Equations | A9f 
|A12| in the absence of the forcing term. Also, instead of using 
an outgoing wave outer boundary condition, we set P = 1 at 
the outer boundary (this corresponds to a choice of normal¬ 
ization). Physically, this scenario would represent a steady 
state oscillation due to irradiation by waves with angular 
frequency oj. The steady state is composed of both ingoing 
and outgoing waves with equal magnitude so that there is 


no net energy transport. It resembles an oscillation “mode” 
of the background structure, although the mode spectrum 
is continuous. 

The solution calculated via this technique is the super¬ 
position of an ingoing and outgoing wave. Each perturbation 
variable can thus be expressed in the form 


Cr(r) = P(r) 


i f,r krdr I —i fn k^dr 

e Jo ^ + e Jo ^ 


= B{r) cos 


rr 

Jo 


krdr 


(B9) 

(BIO) 


where B{r) is a wave amplitude and kr is the radial wave 
number (Equation |14[). In the WKB limit this implies that 
the wave amplitude is 



(Bll) 


which is a smooth (non-oscillatory) function of radius and 
can be calculated from our numerical wave solution. Both 
the ingoing and outgoing energy flux at any point where the 
wave is in the WKB limit is 


E{r) ~ 2puj^r^CaB^ . (B12) 


At any point within the star, we can then define a wave 
leakage timescale 

tleak(r) = , (B13) 

E{r) 


where E{r) is the net wave energy contained below the ra¬ 
dius r 

E(r) . £ f , (BI4) 

and dE/dr given by Equation |B5| 

The quantity tieak(r) is a smoothly varying function 
that describes the amount of time it would take waves to 
leak out of a region below radius r in the absence of an in¬ 
going wave flux. Therefore the value of tieak calculated in this 
toy problem serves as a good proxy for the wave leakage time 
for our forced oscillation calculations. In the absence of wave 
reflection within the star, the value of tieak would simply be 
the wave crossing timescale tcross of Equation |20[ However, 
the wave energy is concentrated within the PNS due to its 
reflecting edge and so in general tieak(r) > tcross(r) because 
the wave energy only gradually leaks out of the PNS. 

To estimate the timescale on which waves damp via GW 
emission, we compute a GW damping timescale 


iGw(r) 


Ejr) ^ 
Ecwir) 


(B15) 


where E{r) is the wave energy contained within radius r, 
and Eowir) is the GW energy emission rate from Equa¬ 
tion The waves become strongly attenuated when their 
GW damping timescale at a radius r is smaller than the 
time it takes a wave to propagate past r. Except for high 
frequency waves (/ > 2 kHz), the GW damping timescale 
is long, tGw(r) > tieak(r). Therefore, GW emission is not 
likely to be the dominant source of wave damping, and only 
a small fraction of the energy contained in the fluid motions 
will be converted to GWs. 
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B3 Wave Amplitudes and Non-linear Effects 

To formally calculate the wave amplitude at a radial loca¬ 
tion r and time t, one must integrate the response per unit 
frequency via Equation |10[ If the wave response is sharply 
peaked at certain values of w, as it is for frequencies near the 
PNS oscillation modes, one can approximate the response 
due to these waves as 

~ (B16) 

where Aa; is the width of the frequency peak in the com¬ 
puted response. This approximation is only valid at times 
tAcu <C 1 and at radii corresponding to fcross(r)Atj ^ 1, 
otherwise waves of different frequencies will deconstructively 
interfere with one another. 

We use Equation [BT^ to estimate the amplitude of 
waves with frequencies near the PNS quadrupolar f-mode 
shortly after bounce, as shown in Figure The value of 
is of order kilometers in the inner 100 km of the super¬ 
nova, translating to displacements of ~ 20% the radius of 
the PNS. 

At these large amplitudes our linear calculations begin 
to break down. We therefore also plot the value of |fcr.^r|. 
Modes are strongly non-linear and are expected to quickly 
dissipate when \kr^r\ ^ 1. At larger radii, the waves become 
increasingly non-linear, which will lead to non-linear wave 
breaking if the waves make it that far. It is also possible 
that the waves are dissipated by neutrino damping or tur¬ 
bulent dissipation before they are able to generate non-linear 
wave breaking. Nonetheless, the fairly strong non-linearity 
(l^r^rl ~ 1) of these waves within the PNS indicates that 
non-linear processes such as three-mode coupling may be im¬ 
portant for these waves. One common outcome of such cou¬ 
pling is the transfer of energy to waves with / ~ 2/f_niode. 
The GW spectra of fully non-linear simulations (K15) con¬ 
tain a peak near twice the f-mode frequency, indicating that 
non-linear effects may be at play. 
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